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THE REFSAT APPROACH TO LOW-COST CPS TERMINALS 


J. W. Sennott, Howard University 
A. K. Choudhury, Howard University 
R. E. Taylor, NASA/GSFC 

ABSTRACT 

This paper describes a concept utilizing a geostationary rererence satellite 
(REFSAT) that broadcasts navigation - aiding signals to low-cost civil-user ter- 
minals which employ the constellation of 24 NAVSTAR Global Positic 15 Sys- 
tem (GPS) satellites for position determination. 

The low-cost receiver design approach eliminates requirements for delay- 
lock tracking and Costas loops; relaxes local oscillator long-term stability require- 
ments from ±0.01 parts per million (ppm) to ±10 ppm; and reduces required 
digital processor capability substantially - all of which lower production costs 
compared to military counterparts. A design goal is to achieve a civil-terminal 
receiver cost of SI 000/set, or less, in volume production. 

This paper describes the signal acquisition, tracking and nosition-fixing 
properties of such a low-cost, dual-channel, L-band, civil user receiver designed 
to receive both GPS navigation and REI'SAT navigation-aiding signals. 
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THE REFSAT APPROACH TO LOW-COST GPS TERMINALS 


I. INTRODUCTION 

Present plans call for world-wide deployment, by mid-1985, of a constellation of 24 NAV- 
STAR Global Positioning System (GPS) satellites. Many excellent articles^ are available detail- 
ing the satellite constellation, as well as the various classes of user equipments under development. 
Although the impetus for GPS has come primarily from the Department of Defense, the signifi- 
cance of a high-accuracy position and velocity determination capability, in three dimensions, has 
not gone unnoticed by the civil sector, specifically the Department of Transportation and the 
National Aeronautics and Space Administration (NASA). Various efforts are underway to ascer- 
tain the potential of GPS for aeronautical, space and surface navigation functions. A theme com- 
mon to these efforts is the development of “low-cost” user equipment concepts which would en- 
able early and wide-spread civil utilization. 

In a recent article by T. A. Stansell,"^ one finds a good summary of the issues involved in 
realization of low-cost civil GPS equipment. As pointed out, one possible avenue to cost reduc- 
tion is alteration of the present GPS signal structure, as realized by a modification to GPS 
spacecraft. 

This paper describes the geostationary reference satellite (REFSAT^ approach having the fol- 
lowing distinct advantages: 

1. No impact on DOD spacecraft or user equipment, while making fun use of existing GPS 
C/A “coarse acquisition” signals. 

2. Low civil-user terminal cost with potential for under $1 000/set given large production 
quantities and Large Scale Integration (LSI) technology. 

.L Provision of acccss/denial command system for civil community. 

The REFSAT approach employs narrow-band carrier-acquisition and frequency-shift-keyed 
(FSK) data signals sent from a geostationary space platlorm, in a configuration depicted in Figure 
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I. An L-Band downlink transmission is sent at the lower edge of the GPS LI (1575.42MHz) 
channel, and is passed by a common receiver front-end into dual-channel IF amplifiers. 

Development of the REFSAT concept involved a careful examination of the signal processing 
functions performed by the conventional NAVSTAR GPS user terminal. A contrast between the 
REFSAT and conventional receiver functions are described in the simplified block diagram. Fig- 
ure 2. 

II. OVERVIEW OF THE GPS-REFSAT TERMINAL CONCEPT 

Before examining the REFSAT user terminal hardware and software requirements it is worth 
while reviewing the GPS signal structure available to civil users, together with the requisite func- 
tions performed in the conventional terminal. 

The operational GPS constellation will consist of 24 satellites in circular 10,900nmi orbits, 
in three 63® (or 55®) inclined planes of 8 each.^ At least 6 satellites are in view at any one time, 
from any point on earth; on the average nine satellites are in view. Each contains an atomic time 
standard which, among other things, generates a periodic 1023 PRN reference C/A code with 
period of 10"^ seconds. A spread-spectrum, L-band, ranging signal is derived by employing bi- 
phase modulation. All satellite clocks are kept synchronized by ground station tracking. Each 
satellite has its own unique PRN sequence, and in this way the two-MHz-wide C/A channel works 
in a code-division multiple access mode. For military users, this has the important advantage of 
simultaneous ranging on two or more signals while offering anti-jamming protection. 

In Table 1 , below, we have summarized the various operations to be performed iii deriving 
a position fix. These are grouped by major function: (1) signal acquisition, (2) signal tracking 
and (3) position fixing. In implementing each, several sub-functions are performed. We ^ ill now 
contrast the GPS and GPS-RFFSAT terminal sub-^fimetions. 
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(a) CONVENTIONAL NAVSTAR GPS USER EQUIPMENT 
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Figure 2. REFSAT vs. Conventional Receiver 











Table 1 

GPS Receiver Functions Comparison 
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In the conventional GPS terminal, signal acquisition begins with onboard computation of 
GPS satellites in view. To accomplish this, the terminal must know, approximately, its own loca- 
tion and have available a recent version of the GPS almanac. For acquisition to go forward, Dop- 
pler estimates must then be computed for the subset of GPS satellites having good geometry. 

This leaves the GPS demodulator, typically, with a-priori GPS carrier frequency uncertainty of 
several hundred hertz (Hz). However, to stabilize the receiver front-end local oscillator to this 
accuracy requires an oven-stabilized L-band synthesizer. To complete the signal aojuisition func- 
tion, the GPS receiver generates a local C/A code reference for each selected GPS satellite, and 
searches for the cross-correlation peak in 0.5-1 .0 microsecond steps, spanning the 10*^ second 
code period. 

The acquisition function, as performed in the GPS-RhFSAT user terminal, differs from the 
above in several important aspects. Before examining these differences, the REFSAT down-link 
signal format. Figure 3, must be understood. The REFSAT signal structure consists of an FSK 
data stream divided into 128 bit sub-frames. Each contains a Doppler gradient word, C/A code 
select word and satellite position (referenced to REFSAT system frame epoch). The overall FSK 
transmission is centered on frequency ccrj. Considering the GPS-REFSAT terminal acquisition 
function, the aided terminal first acquires the FSK data stream, which is accomplished with a 
non-coherent FSK demodulator. Since ihe REFSAT concept employs a dual-channel L-band 
receiver with a common L.O., drift effects are common to both the GPS and aiding signal. Thus, 
a long-term stability as poor as ±10ppm may be employed with no effect on Doppler acquisition. 
Also noteworthy is the satellite selection procedure which, in the GPS-REFSAT terminal, is de- 
veloped directly From the FSK data words. The GPS-REFSAT terminal completes the signal 
acquisition function with a C/A code scan identical to thai ibr the conventional terminal. 

Consider now the .dgiial tracking (unctions listerl m Table I. In (he conventional GI’S ter- 
minal signal tracking occurs m three stages. Firstly, the code delay must be (racked .nid an 
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Figure 3. REFSAT Data Format 



aligned or “punctual” code stream generated. This is impIemt.Tted \vith either a delay-lock loop 
or tau-dither technique. Regardless, the end result must be a continceasly tracked code genera- 
tor with delay error held to within one-tenth of a microsecond or better. Secondly, Doppler un- 
certainty must be eliminated. This may be accomplished by stepping the frequency synthesizer 
and measuring correlator output. Once the Doppler uncertainty is reduced to 10-20Hz, carrier- 
phase and GPS data is recovered using a Costas loop and the aforementioned punctual code. 

The GPS-REFSAT terminal performs the signal tracking function in a substantially different 
fashion. A continuously-tracked punctual code is not necessary. For most users a code-delay 
interpolation scheme is sufficient. With this technique, described in Section 111, code step pairs 
separated by one chip are employed. This simplifies the code generators and tracking filters. 
Further, the Costas carrier phase tracking loop is eliminated. 

The user terminal function performed last is position fixing. As shown in Tabie 1 . this is 
composed of three sub-functions. In the conventional terminal ephemeris parameters, contained 
witliin the 50 B/S data message, are jodated hourly. The GPS termin<.! employs these parameters, 
together with a software model for the satellite orbits, to compute satellite coordinates in real 
time. This requires substantial memory and compui'.lional resources. Given time delay (and Dop- 
pler) estimates, the GPS receiver then computes portion, time, velocity and bias variables of 
interest. 

The GPS-REFSAT terminal, on the other hand, offers a simplification. The REFS AT FSK 
data stream (Figure 3) provides, every four seconds, updated GPS satellite coordinates. This 
eliminates most of the storage and computation associated with hourly updating of GPS satellite 
ephemeris parameters. In Section V, we will discuss the position fix and satellite co-cordinate 
update approach in more detail 

At this point, it is useful to present simplifieo block diagrams for the con.'.ntional and aided- 
user terminals. Figures 4 and 5 rcspectiveh For the conventional terminal we have assumed a 
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Figure 4. Simplified Block Diagram Un-aided GPS Terminal 















‘*single-chaiuier’ version in which acquisition and tracking circuitry is shared, sequentially, be- 
tween different GPS satellites. Components in the GPS receiver include RF/IF hardware, L- 
band synthesizer and frequency reference, code tracking loop, carrier tracking loop and digital 
processor. For the GPS-REFSAT receiver the L-band synthesizer and reference have been elimi- 
nated, along with punctual code and Costas loop components. Further, the code generator and 
digital processor have been simplified. 

In the following sections a more detailed description (and analysis) of the acquisition, track- 
ing and position fixing functions is provided. 

III. SIGNAL ACQUISITION 

At the user terminal a given GPS satellite signal will be buried in both additive Gaussian 
noise and signals from six to eight other GPS satellites. The desired signal is received with un- 
known Doppler and range delay. Ultimately, the receiver must compute the relative Doppler and 
delay between this satellite and all others which may be visible. The initial acquisition function 
is carried out with an M-ary hypothesis test, over the appropriate range ti.Tie-aeIay and Doppler 
f«eqi»ncy cells. In mechanizing the above test it becomes nccess;iry to generate, locally, time 
arid frequency-shifted C/A codes. One then generates the in-phaNe and quadrature correlations, 

I and Q. Under microprocessor control, a family of statistics are generated by sweeping out the 
total Doppler and delay uncertainty. These non-coherent statistics are denoted by A(i J). where 
i and j are code and VCXO indexes, respectively. 

Under noiseless conditions, a given A(ij) appears in Figure 6. A maximum corr ition value 
results when the received delay and Doppler parameters match those generated by the local code 
and voltage controlled crystal oscillator (VCXO). Correlation on the time delay axis is essentially 
zero beyond 10"^ seconds offset, with the Doppler zero-crossing at 10^ Hz. With channel noise 
present a delay-Dopplcr pair (ij) far removed from the peak ma> yield maximum correlation. 
The probability of such an event, computed in the Appendix, is denoted as |li|. To cembat 
this type of error, a sequential test is employed. 
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This approach 'jses a 'Tun ^ ^mparison” threshold. Statistic A(i,j) is compared with A(i-lJ), m 
times. If A(iJ) > A.(i-I,j) in each test, a “hit” is registered. The process is repeated until the 
total set of dela}’ Lfidexes is searched, or a “hit” is registered. If the number of delay steps to be 
scanned is s, th?i' he probabilities of interest are computed to be; 


Pj (start false tracki = s (1/2)™ 

(0 

Pj [detect true delay] = (1 - P||(E)]™ 

(2) 

Pj [dismiss true delay] = 1 - [1 - P||(E)]™ 

(3) 


Doppler uncertain affecting values of P||(H), expressed by the Marcum-Q function defined in 
the Appendix is a predominant factor. Figure 6 shows the rapkl fail-off of correlation with in- 
creasing frequency offset. The entire delay index scan is repeated for each 200 Hz - increment 
of Doppler uncertainty. 

The obser ed Doppler will be different at each receiver. For the CONUS, total Doppler 
span for each satellite is approximately 4kHz. Typical Doppler contours are shown in Figure 7; 
the gradient is ?*jout one Hz/km. Now if the receiver has some a-priori knowledge of position, 
say to ilSOkm, an .nitial predict '.>1 of Doppler to several hundred Hz can be .made. However, 
a single-channel receiver’s L-band stability must be this accurate as well since a prolonged Dop- 
pler scan would be necessary. As discussed previously, the R6FSAT dual-channel receiver solves 
this problem. Also, the REFSAT carrier can be transmitted with a known frequency offset rela- 
tive to the GPS satellites of interest. The received frequency offset may then be computed by 
the user, by making use of th«^ Ooppler gradient word and an initial position estimate. The fre- 
quency detuning error (Indian per second) as seen at the receiver is: 

- Xi * 1 "1 r ’ ^Rs 1 

" IaI — 0)LOJ ■ Vlo "IrsI ' <^1.0 

liq. (4) reflects the R17iF struclurc of Figure 5. As previously noted, the term cancels be- 
cause of the dual-cliaiii I .cceiver. Additionally, for many users, the user-lo-REFSA T range-rale 


1 .^ 
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term, |Fq - F|^|, is small. The initial value of j is then found from; 

jAu) = |ro ‘-I aI (5) 

C 

where the term on the right of Eq. (5) may be estimated from the REFSAT Doppler gradient 
word for GPS satellite A. 

Returning now to Eqs. (1-3), let us assume that the Doppler detuning error has been reduced 
to within 150 Hz. Considering delay-detuning, the worst-case condition is clearly for an arrival 
time which is exactly between two delay steps, and on the average the error is 0.2SAr. Now go- 
ing to the 43dB-Hz performance curves. Figure 1 1 of the Appendix, the false hypothesis prob- 
ability Ph(E) is 0.008. Using a run repeat index value m = 10 yields a “start false track” prob- 
ability of 0.976 and a “detect true delay” probability of 0.93. This is for the worst-case in 
which all 1023 delays must be scanned before getting to the correct delay. In this situation one 
false track may be passed to the track verification phase but, with probability 0.93, the true 
track will be included as well. The overall scanning process, for m = 10, requires about 20 sec- 
onds per satellite. Total acquisition time for 4 GPS satellites would be on the order of 80 
seconds. 

A hybrid analytical/simulation model will be developed to obtain more exact performance 
statistics involving received signal strength and the parameters s and m. 

IV. TRACKING AND INTERPOLATION 

Unlike the acquisition mode, wherein many A values may have to be scanned, once tracking 
commences index gates are continually available for each of the satellites A, B, C, and D. These 
are maintained in a repeat-scan fashion whereby codes A, B, C, and D are checked sequentially 
using the two closest values of step index, say i and i + 1. In ail. eight statistics are generated 
and checked during each cycle. The process is repeated every 1 2 milliseconds, depicted in Figure 
8. An interpolation scheme is then employed to generate the time delays used in position 
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computation. As pointed out in Section II, this can be done since a punctual code is not re- 
quired. For satellite A, the time delay estimate is given by; 

"A='*AAr + lA (6) 

AaUJ) > AaO + l,j); 

The interpolation function I^ is composed of the difference of two adjacent test statistics, nor- 
malized by a gain constant formed from their sum. In Figure 9, we have plotted I^ and A, for 
noiseless conditions, as a function of code delay offset Ar^. Iiuci eolation standard '‘‘Aviation is 
developed in ihe Appendix and tabulated below for a range of carrier-to-noise (CNR) ratios: 


CNR 

At = lO"* sec. 

At = 0.5 X 10-^ sec. 

(40dB-Hz) 

120 m. 

60m. 

(43dB-Hz) 

75 m. 

37 m. 

(46dB-Hz) 

60 m. 

1 

30 m. 


It is significant to point out that the above “raw” errors can be reduced substantially, for 
most civil-user applications, by post-interpolation smoothing. Furthermore, the observations can 
be updated rapidly, at a maximum rate of about 80 times/sec., making them suitable for certain 
high-dynamics user applications. 


IAa(U)J‘'^ -IAA(i+ At 

IAa(U)1'/2 + [Aa(1+ 1J)I'^J 2 


where, for 



A typical signal processing sequence, controlled by the microprocessor, is depicted in Figure 
8 for four, visible, GPS satellites. Implicit in this sequence are the following operations; 

• Recent interpolations are averaged. 

• Interpolations far removed from average are edited out. 

• Code phase step is incremented or decremented when the average interpolation is outside 
sivcil'icd limils. 
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(a) CODE-DELAY INTERPOLATION FUNCTION 
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(b) CORRELATION-INTERPOLATION COMPONENTS 
Figure 9. Code-Delay Interpolation F'unction and Correlation-Interpolation Components 
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V. POSITION FIXING 


The REFSAT-aided receiver must, as does the conventional GPS receiver, convert estimated 
code arrival times into user n^jiuon. In e*»her «^se, the solution depends upon a knowledge of 
GPS satellite positions at the instant of ranging code transmission. The GPS-REFSAT receiver 
microprocessor keeps track of these positions by decoding the RFFSAT FSK data stream (Figure 
3), whidi contains direct coordinate information referenced to the GPS C/A code period emitted 
at event time ^RS- The REFSAT approach eliminates the need tor satellite dynamical models, 
and associated computations, normally required. The method for obtaining a user position fix 
will now be discussed in more detail. 

In visualizing the various event times, it is helpful to view each C/A code period of 1023 
chips as a single pulse, emitted by the GPS satellite at the onset of that period. We neglect in 
this discussion the small clock offsets which exist between different GPS satellite clocks, and as- 
sume that all codes are synchronized. Every four seconds the REFSAT system defines one of 
these synchronous emission times as having occurred at T^rs- If but one set of pulses were emit- 
ted by the GPS satellites every four seconds, and if the receiver could measure the arrival times 


on parallel channels, the signal arrival times for GPS satellites A, B, C and D would be: 

’•a =IIo(^a)-Ia(Trs)I/C + Trs (8) 

*"8 ~ llo ~Ib^T^Rs)I/C ^RS 

~ llo ■*" ^RS (10) 

Td = llo (^d) -Id(Trs)I/C + Trs (11) 

The observed delay estimates,"?), are equated to the above theoretical delay values to get: 

^A =llo(rA)-'-A(TRs)l/C + TRs (12) 

~ llo ~ (Trs)I/^ Trs (13) 

“ llo (^(') ~Ic(Trs)I/(" Trs (14) 

" llo (^d) " ‘'d(Trs)I/(^ )^RS (1^) 
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In solving these the (usual) assumption is made that fg, user position, is constant between pulse 
arrival times, and that the GPS satellite positions at the reference time are accurately known. The 
latter information is encoded directly with REFSAT data blocks A-D (Figure 3). Each coordi- 
nate is encoded with 23-bit accuracy, for an error of about 2 meters. In solving for^Q. the 
microprocessor also provides T{^. A time-of-day word encoded at the start of each data block 
may then be used to set the user clock with that of the ground stations. 

The idealized scheme above is complicated by the fact that GPS code periods are emitted 
every millisecond, with receiver interpolation accomplished sequentially over a 12 millisecond in- 
terval. Further, the user may want position fixes anytime during each frame. In such cases, the 
receiver software must account for GPS satellite motion from positions at T{^. This is a straight- 
forward procedure, assuming the user position is known to within several hundred km. The “off- 
set” satellite positions, at k| code periods following T^^, are estimated to be: 

ti (Trs + kjAT) = li (Trs) + kiATti (16) 

Satellite velocity estimates?} of sufficient accuracy result simply by differencing satellite position 
between data frames. Worst-case satellite positioning error is on the order of 7 meters, plus any 
uncertainty due to ground control station tracking. 

The major position-fixing error sources have been tabulated in Table 2, for three different 
position fix rates. For fixes smoothed over one second the delay estimation contribution is 6.7 
meters, whereas for \ery high update rates this term contributes as much as 60 meters. For most 
civil users, the dominant error source will be due to ionosphere propagation. A 35-meter value 
is worst-case, and about an order of magnitude larger than experienced at most times. Neverthe- 
less, it may be desirab'e to include a propagation correction word in each of the REFSAT data 
blocks. Analogous to the Doppler gradient word, this would allow a direct correction in real- 
time, given approximate knowdedge of user position. Siicli a .scheme should allow porition I'ixini; 
to substantially better than 100 meters. 



TabJe 2 

Erroi Budget Summary 


Error Terms 

fix rate^ * ^ 
I per sec 

fix rate^^^ 
10 per sec 

fix rate 
80 per sec 

(2) delay estimation 
(40dB-Hz CNR) 

6.7 m 

21m 

60 m 

(3) GPS satellite 
positioning 

9 m 

9m 

9m 

(4) propagation 
uncertainty 

35 m 

35 m 

35 m 

(5) GDOP multipUer 

3:1 

3:1 

3:1 

Total fix error 

152 m 

195 m 

312m 


Notes: 


(1) assumes constant velocity and independent range-rate solution, 

(2) interpolation perfoimed with 1/2 chip offset, 

(3) 2-m error contribution from GPS Master Control Center and 7-m from REFSAT FSK encoding, 

(4) may be reduced with REFSAT propagation correction transmission, 

(5) Geometric Dilution of Precision (GDOP). 


VI. FUTURE ACTIVITIES 

NASA is planning to demonstrr***, operationally, the REFSAT concept during CY 1979; a 
“breadboard” REFSAT user-terminal receiver is currently being designed for this purpose. 
NAVSTAR GPS Phase I satellites (6 expected in orbit by Fall-1979) will be employed to deter- 
mine geographical location with the breadboard REFSAT receiver. 

The operational demonstration of the REFSAT concept will be backed up with a software 
simulation and analysis model for determination of user terminal position error. 

By mid--1979, a detailed cost analysis of a REFSAT user terminal design will be completed 
to ascertain if low-cost production cost goals can be achieved. 



VII. CONCLUSION 


The REFSAT-aided GPS concept appears to offer significant advantages tl tribute to 
low user terminal production cost, as evidenced by reduced electronic circuit complexity, and 
relaxed-component precision. 

The REFSAT approach eliminates delay-lock tracking and Costas loops, relaxes receiver 
front-end local osdllator long-term stability requirements; and reduces microprocesi'or require- 
ments by eliminating satellite almanac and satellite positioning computations-all of which lower 
production costs compared to militap' counterparts. The GPS C/A code navigation performance 
capability has not been sacrifu-.d, making the navigation capability, for example, more accurate 
than LORAN-C or VORTAC. 
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APPENDIX 


TEST STATISTIC ANALYSIS 


The REFSAT user terminal genarates mm-«oherent opde conelation test statistics makuDQ 
use of the signal piocessing anpioach discussed previously. Figure 8. Such statistics are employed 
both during Mgnal acquisition and interpolation (tracfcii^). The Ibliowing analyas is similar to 
that for non-coheient radar and non-coherent orthogonal data systems. 

Cbnsider first an equivalent noise model for Figure S. as Figure 10. All GPS channel noto, 

* 

including that from tlie front-end .-band mixer and W'XO mixer, is replaced by a single noise 
source of Nq watts/Hz. The noise level is referenced, to a. unit -amplitude carrier. Tlie Doppter 
and delay detuning errors between the two ^gital i^aths are and = r - t\ respec- 

tively. In-phase (I) and quadrature (Q) channel outputs followiiig integration arc: 


/ V ^ ^ r /- 

•^c<»s|u)t + lit - T> ♦ dl + W,.tt)| ^«»s|«j't + PtI - t') -t PI dl {I7» 

*■ 


rtj+AT , 


*-s ~ 


cosicat + P(t - r) + pi + W, 


•^sinlcj'i 


t + P(t -T) + Pl|^dt (18) 


In forming the sum-square statistic from these quantities, the phase error term p - p' is eliminated. 


We first observe that both and are Gaussian random variables. Working on their ex- 
pectations, and replacing tlie 0, rr function !*(•) with a function P'(*) which takes on values tl, 
we have: 


l-.ll., ItJj., Ti.| 


FlUilco^. T^.| 



tj + AT 


tj +AT 


P'(l -t)I**(I -t'P'os|u>1 + p|cos|tA)'l + p'|dt 


P'(t - t)I*'(I - t'K'c ; + plsiii|cj'{ + p’Idl 


C 


(20) 
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We defiiie function R(*), the autooonehtk>n of the seqiwnce ?*(•) as. 


R(t-t*) 



+ NAr 

P'{t-T)P'(t-r')dt, 


( 21 ) 


whioe N is the PRN code period and Ar is the duration of one “chip”. The product function 
P'(t -r)»P'(t - r') contains spectral components at multiples of the baric PRN period. In the 
case under sAuiy these components riart at DC and run to several MHz. Their spacing is approx- 
imately one kHz. The non-^X? terms are ctasrifled as “self-^ioise.” It is easy to show that the 
DC term is given by (21). For large processing gain (i.e., integration periods on the order of one 
code period or larg») we may neglect all but this term. One then has for the above expectations: 


EILsIcOe.Te) = RIt^) 


E(Lclc<>e,r^] = R(rg) 



+ AT 

cos[(ot + 0]cos[(c't + 0'ldt 

+ AT 

cos[<ct + 0]sin[co't + ^')<lt 


(22) 


(23) 


We express the integrands in terms of difference and double frequency terms. For AT = 10*^ sec 
and (o-2n X 10’ rad/sec, the double frequency term may be neglected. After integration we 
have; 


R(fe) ( 'I 

E[Lsl<OQ,re] = 1/2 -Isinla + WgAT) -rin«> 

We V J 

R(’’e) r ^ 

ElL^>|We,Te' • 1/2 •< cost® + WgATl - cos® > 

w- > ) 


where 


o = Wgti + (0 - 0') 

In a similar fashion we may compute the (equal) variances for and as, 


Nc 

2 


(24) 

(25) 


( 26 ) 
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We can show tint the two random vanaUes are statistically independent by slnming orthogomditv 
of the in-phase and quadrature reference signds in the interval {t|, tj AT] ; 


rt+AT 


A 


OOS(w't + P(t - t') + d']^[(o't + P(t - t') + p’ldt SS 


(27) 


+ AT 


1/2 I sin[2w’t + 2P(t-T') + 2d']dt«0. 

The do»ity function for the (root) sum-equare test stati^ic is then readily shown to be^; 


P^%( 2 ) = ^ exp - {/.2 + ij2l/2<r2 lo 


(s) 


(28) 


where 


ij s ^ElLsIWg.Tel* + Etl^lWe.Tel^ 

Now returning to the acquisition problem, the basic event probability of interest is: 

P„[E] = P,[A(iJ|Te,cOe) > A(i*J|Te, We)l 

where i* is the code delay step which minimizes r^, and i is a code delay step earlier or later 
than i*. This probability can be developed using the density of F.q. (28), or by a derivation 
from Van Trees^. Using the latter we have: 

PIxf + x| > x2 + x|) = 1/2H - Q0,a) + CKot,/5)l 

where (X,) is the Marcum-Q function, Xj are Gaussian, and 

1 1/2 


(29) 


(30) 




Applyini: this dirci'lly to l-q. (2'1), wo nolo t'roin i-'q. (24) ami tlio sharp autocorrelation function 
R(*), that a = 0. Further, combining Eqs. (24). (25) and (26) we have for 


26 



( 31 ) 





2 


Ng 


AT 


Apidyii^ standard ideniitks for the Marcum-Q function we then have, 

Pu(El = U2e-P^ (32) 

ThB fimctk>n k shown in Figure 11 for a range of delay detuning values, parameterized in CNR 
and frequency d<^£unmg. 


Last^, we wish to discms in more detail the interpolation function of Section IV. The 

quantity desired is the variance of Oi^. For convenience we reproduced Eq. (7) below: 

A 

I . L - - IA«(i ♦ I 

* " \ ' lAAOjHW + IA,(i * lj)lW2J 2 

This is a function of four random variables with densities of the form of Eq. (28). The two 
variables in the denominator, under noiseless conditions and with unit-amplitude received carrier, 
sum to unity as seen in Figure 9b. For the range of (TNR and integration times of interest, it 
is reasonable to assume that the noise contribution in the denominator is much smaller than 
unity. Now in the numerator, we have the difference of two (positive) random variables. For 
= Ar/2 their expectations cancel. The variance oi interest is then, approximately; 



Wc can over-bound Oy^^i/2 (26). 
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Table 3 

Summary of Notation 


AaOJ) 

1" Ua> *b» *c» J 
At 

Ia» ^b» 4:» 
i= Ua>Jb»Jc*JdJ 
A(*> 

WRS 

‘^LO 

We 

Trs 

AT 


Ia 


Ia 


lo 


So 


P(t) 

wc(t) 

Ng 

Lc> Ls 


N 


Test statistic for GPS satellite A generated with delay step i and VXCO 
step j. 

Maximum-likelihood code delay steps for GPS satellites A, B, C, and D. 
Code Delay Off-set Error. 

Delay time step size. 

Delay interpolation constants. 

Maximum-likelihood frequency steps for GPS satellites A, B, C, and D. 
VCXO angular frequency step size. 

Carrier frequency emitted at GPS satellite. 

Carrier frequency emitted at REFSAT. 

Receiver LO frequency. 

Frequency de-tuning error. 

Start time oi REFSAT frame sync word at REFSAT satellite. 

Statistic integration time. 

Positions of GPS satellite A. 

GPS satellite A velocity. 

User Terminal position estimate. 

User Terminal velocity estimate. 

C/A code waveform. 

Noise waveform at GPS Channel IF output, unity amplitude signal 
reference. 

Noise spectral density at GPS Channel IF output. 

In-phase (cosine) and quadrature (sine) integrator outputs. 

Phase error. 

Variance of integrator outputs. 

Number of chips in C/A code. 



Table 3 (Continued) 


R( ) 

Autocorrelation of C/A code. 

a./J 

Error exponenti. 

PhIEI 

False index error probability. 

Qi , ) 

Marcum-Q function. 

0 } 

Delay interpolation error variance. 

c 

Propagation velocity. 
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